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I. INTRODUCTION 

The problem of temporal and spacial variation of the 
fundamental physical constants is actively discussed in 
the literature during last several years. A recent review 
of its current status can be found elsewhere [J]. One of 
the reasons stimulating this activity was the discovery of 
acceleration of the universe (for a review see jll), what 
is usually associated with the existence of the dark en- 
ergy. Modern theories describing cosmological evolution 
predict that the dark energy may cause variations of the 
fundamental constants. 

A statement that the fine-structure constant a has 
possibly changed during evolution of the universe was 
made by Australian group in Ref. (SL Other astrophysi- 
cal groups do not confirm this result 3, , and hence new 
laboratory and astrophysical investigations are required. 

Laboratory studies of hypothetical variation of the 
fine-structure constant are based on the fact that tran- 
sition frequencies in atoms depend on aZ, where Z is 
atomic number. Supposing that the nowadays value of a 
differs from its value in the earlier universe we can study 
space-time variation of a by comparing atomic frequen- 
cies for distant objects in the universe with their labo- 
ratory values. In practice, we need to find relativistic 
frequencies shifts, determined by so-called q factors, ac- 
cording to 



wiab + qx, X = (a/aiab)^ 



1 



(1) 



Most advantageous for these studies are the atoms and 
ions for which q factors of transitions between certain 
states significantly differ from each other. At the same 
time these elements should be abundant in the universe 
to provide sufficient observable data. Fe i and Fe ii dis- 
cussed in this paper meet both these conditions. Spec- 
tra of Fe II were used by several groups [1, 0, Q within 
"many-multiplet method" and by Levshakov et al. within 
"single ion differential alpha measurement method" Q. 
The atoms of Fe i were observed in resonance ultravio- 
let lines in two damped Lya systems at z=0.452 Q and 
2=1.15 [8]. According to @ the spectra of Fe i are also 



observed for the high redshift quasars and may be used 
for the a-variation search. In this work we calculate the 
q factors for the transitions in Fe i from the even- and 
odd-parity states to the ground state. 

As was pointed out in many papers including [1, 
[Til , [l^ . [isj one of the problems that occurs in the study of 
possible a variation is a necessity to separate this effect 
from the isotope shift (IS) effect. A method to resolve 
this problem was suggested in [l2| • This method requires 
(along with a knowledge of atomic relativistic coefficients 
q) precise calculation of the isotope shift coefficients. 

In this paper we carry out isotope shift calculations in 
Fe I and Fe ii. We compare the results obtained for Fe ii 
with other available data. To the best of our knowledge, 
the calculation of the IS in Fe i is performed for the first 
time. 

The paper is organized as follows. Sec. [Til is devoted to 
the method of calculation of the properties of Fe i. We 
discuss the results obtained for the q factors. In Sec. Hill 
the method of calculation of isotope shift is described. 
We present the results of the isotope shift calculation 
in Fe i and Fe ii. Finally, Sec. IIVI contains concluding 
remarks. Atomic units {h = |e| = rrie = 1) are used 
throughout the paper. 



II. CALCULATION OF q FACTORS FOR Fe I 
A. Method of calculation 

To find q factors we need to solve the atomic relativis- 
tic eigenvalue problem for different values of a or, re- 
spectively, for different values of x from Eq. ([T]). We can 
calculate atomic frequencies w± for two values x = ±1/8 
of the parameter x. Our experience shows that such a 
choice of x allows us to meet two conditions. The value 
I a; I = 1/8 is usually sufficiently small to neglect nonlin- 
ear corrections and sufficiently large to make calculations 
numerically stable. The corresponding q factor is given 
by 



q = 4(ll'_|_ — u!-). 



(2) 
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The ground state configuration of Fe i is (Is^... 3c?^4s^). 
Since it fias eight electrons in open shells its spectrum is 
rather dense and complicated. Due to proximity of the 
levels with the same total angular momentum (especially 
in the middle of spectrum which astrophysically is most 
interesting) they strongly interact with each other. All 
this makes calculations of Fe i very difficult. To the best 
of our knowledge the only calculation of q factors for neu- 
tral iron was carried out recently by Dzuba and Flam- 
baum in Ref . [31 . They used a simple method combining 
ab initio Hartree-Fock and configuration interaction (CI) 
technique with some semiempirical fitting of energy lev- 
els. 

In this paper we make pure ab initio calculations in the 
frame of the eight-electron CI method. The [ls^...3j»^] 
electrons are treated as core electrons while 3d, 4s, and 
4p electrons are in the valence space. The number of con- 
figurations accounted for in our calculations is noticeably 
greater than in As a result an effect of configura- 

tion interaction is treated more accurately. Below we will 
discuss it more detailed. 

We started from solving the Dirac-Hartree-Fock equa- 
tions. The self-consistency procedure was done for 
(Is^... 3c?^4s^) configuration. After that the 4pi/2 and 
4p3/2 orbitals were constructed as follows. All electrons 
were frozen; one electron from the 4s shell was moved to 
the 3d shell and another electron from the 4s shell was 
moved to Ap shell. Thus, the valence orbitals 4pi/2 and 
4p3/2 were constructed for the idJAp configuration. 

On the next stage we constructed virtual orbitals. We 
used the method described in [l^ [l^ and applied by us 
for calculating different properties of Fe ii [17|. In this 
method an upper component of virtual orbitals is formed 
from the previous orbital of the same symmetry by mul- 
tiplication by some smooth function of radial variable r. 
The lower component is then formed using kinetic bal- 
ance condition. 

Our basis sets included s,p,d, and / orbitals with prin- 
ciple quantum number n < N. We designate them as 
[Nspdf]. We have carried out the calculations of en- 
ergy levels, g and q factors for three basis sets with 
N = 4-^6. Configuration space was formed by single and 
double (SD) excitations from the configurations 3(fAs^, 
3d^4s 4p, and 3d^4p. 

Additionally we studied Breit corrections including the 
Breit interaction into consideration. For Aspdf basis set 
we computed q factors in the Coulomb-Breit approxima- 
tion and compared them with the results obtained in pure 
Coulomb approximation. As analysis showed for a ma- 
jority of states the Breit interaction changed the values 
of the q factors only at the level of few percent. For this 
reason all results which we discuss below are obtained in 
the pure Coulomb approximation. 



TABLE L Energy levels (cm^^), g factors, and q factors 
(cm"'^) for the (Sd^is^) a^Dj states. The values are obtained 
for the [Gspdf] basis set. 
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B. Ground multiplet ^Dj 

The ground state FS transitions in mid- and far- 
infrared were observed in emission for different redshifts 
for a number of atoms and ions (see,e.g.,[ll, [H, |20|). 
The infrared FS lines of the neutral iron have not been 
detected yet in astronomical objects. But a such detec- 
tion is expected in extragalactic objects at a new gener- 
ation of telescopes like the Stratospheric Observatory for 
Infrared Astronomy and Far Infrared and Submillimeter 
Telescope. 

We start discussion from the results obtained for tran- 
sitions between the states of the ground multiplet. The 
CI space corresponds to SD excitations from the config- 
uration 3d^4s^. In Table U we present results obtained 
for the [6spdf] basis set for energy levels, g factors, and 
q factors for transitions from the ground state ^D^. 

A method to use these fine-structure (FS) transitions 
to study a variation at very high redshifts was suggested 
in Ref. 22]. This method crucially depends on the differ- 
ences of the dimensionless sensitivity coefficients defined 
as 

Qj,j-i = {qj - qj-i)/^J,J-i, (3) 

where ujj^j^i = Ej — -Ej-i is the frequency of the FS 
transition J ^ J — 1. 

It was shown in [22] that for the levels of the ^^~^^Lj 
multiplet the difference, AQ, between the dimensionless 
sensitivity coefficients and is given by 

the following formula 

AQ ^ Qj,j-i-Qj-i,j-2 - ^ ( ) -1, (4) 

which links AQ to the experimentally observed FS tran- 
sition frequencies. This analytical expression is valid up 
to the terms of the order of (aZ)^. Note, that in the first 
order in the spin-orbit interaction the left hand side of 
Eq. (jll turns to zero. For this reason AQ is very sensitive 
to the ratio of the FS transition frequencies. 

As is seen from Table U we reproduce the FS intervals 
with a few percent accuracy. The uncertainties of the q 
factors listed in Table [Jean also be estimated at the level 
of a few percent. Such inaccuracies in calculation of the 
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TABLE II: AQ for the FS lines within the ground multiplet 
a'^Dj. Numerical results are calculated using Eq. Q for the 
basis set [6spd/]. Analytical expression Q is applied to the 
experimental frequencies as in Ref. j22| . 
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FS transition frequencies and q factors, though small, are 
sufficient to result in significant differences between the 
values of AQ calculated with use of Eq. ([3]) and obtained 
from Eq. (g]). 

To illustrate it we present in Table |TT] the differences 
of the sensitivity coefhcients (AQ) of the FS lines within 
the ground multiplet a^Dj obtained by both methods 
mentioned above. The corresponding entries in Table |TT] 
are denoted as Eq. ^ and Eq. In the first case wc 
used the q factors and the calculated transition frequen- 
cies from Table H] In the second case the Eq. ^ and the 
experimental frequencies were used. 

Comparing these results we see, that there is only qual- 
itative agreement between numerical and analytical ap- 
proaches. For instance, Eq. @ predicts decreasing of AQ 
with decreasing J. When we use Eq. ([3]) we observe the 
same behavior. At the same time, quantitatively, a dis- 
agreement is rather significant. We think that numerical 
errors in computing the q factors and the FS transition 
frequencies are more essential for calculation of AQ than 
the terms ~ {ctZY (and higher) neglected in Eq. For 
this reason we consider the values obtained from Eq. @ 
as more correct. 



C. Odd parity levels 

The transitions from the ground state to the odd- 
parity states are observed in absorbtion in the spectra 
of quasars. The results of the eight-electron CI calcula- 
tions of the energy levels, g factors and q factors for the 
odd-parity states of the 2><PAsAp configuration of Fe i are 
listed in Table IIIII We have carried out calculations for 
the three basis sets [(4 Q)spdf]. The CI space for each 
basis set corresponds to SD excitations from two con- 
figurations idf'AsAp and idJ Ap. The agreement between 
theoretical and experimental energy levels was system- 
atically improved with increasing the basis set and best 
results were obtained for the longest \6spdf\ basis set. 
In Table|TTT]we present the results obtained for the largest 
[Qspdf] basis set. We restrict ourselves by consideration 
of the states with total angular momenta J=3,4, and, 
5 because only to these states there are strong electric 
dipole transitions from the ground (3(i^4s^ ^I?4) state. 

Comparing the q factors found in this work with the 
results obtained in Ref. we see reasonable agreement 
between them for majority of the states. At the same 



TABLE III: Energy levels (cm~^), g factors, and q factors 
(cm"'^) for the odd-parity energy levels of the ZdfiAsAp con- 
figuration. The values are obtained for the \Qspdf] basis set. 
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time some differences reach 40%. One of the reason of 
these discrepancies is the strong configuration interac- 
tion for ceratin levels that significantly influences on the 
q factors. For instance for the astrophysically interesting 
X and x ^F^ states of the 3(i^4s4p configuration we 
obtained in one-configurational approximation the fol- 
lowing q factors: q{x^D^) = 835 cm~^ and q{x^F§) = 
861 cm^^. As is seen from Tabic Hill these values are 
more than two times smaller than those obtained for the 
[Gspdf] basis set, when a large number of configurations 
was included into consideration. At the same time we see 
that even for the largest basis set the number of config- 
urations taken into account is still not sufficient to cor- 
rectly reproduce the order of certain states. For instance, 
the energy levels of the x ^Dj multiplet lays higher than 
the energy levels of the x^Fj multiplet. Since this range 
of spectrum is very dense and there are many nearby lev- 
els with the same total angular momenta it can lead to 
incorrect interaction of the mentioned states with their 
neighbours. 

Now we will discuss the odd-parity states of the Sd'^Ap 
configuration. Correct calculation of different properties 
of these states is more difficult than calculation of the 
states belonging to the 3c?^ 4s4p configuration. First, the 
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TABLE IV: Energy levels (cm^^), g factors, and q factors 
(cm~^) for the states of the 3(fAp configuration. The val- 
ues are obtained for the 5spdf basis set. Results obtained in 
Ref. are given for comparison 
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states of the 3d^4p configuration are located higher in en- 
ergy than the majority of the states of the 3(fi AsAp con- 
figuration considered by us. Second, the energy levels of 
the 3(fAp configuration belonging to different multiplets 
are located very close to each other and, respectively, 
strongly interact to each other. An additional problem 
is that the HFD equations were solved self-consistently 
for the configuration SdPAs^. The configuration 3(f Ap 
differs more significantly from it than the configuration 
3dP 4s4p, in particular, by the presence of an extra d elec- 
tron on the 3d shell. As a consequence, it is more difficult 
to reproduce the correct transition frequencies from the 
states of the 3d^ 4p configuration to the ground state. As 
our analysis shows these calculated frequencies tend to 
be larger than the experimental ones. 

To investigate how different properties of these states 
will change when the basis set is increased from [4spdf] 
to [6spdf] we have carried out the calculations for all 
these basis sets. The agreement between theoretical and 
experimental frequencies was at the level of 5-10%. As 
it turned out the best agreement was achieved for the 
[5spdf] basis set. This obviously indicates that for this 
configuration we are still far from saturating the CI space. 
In Table IIVI we present energy levels, g and q factors 
obtained for this basis set. 

As is seen from Table IIVI the agreement between the- 
oretical and experimental energy levels is worse than it 
was for the states of the Sd^isAp configuration but nev- 
ertheless the largest difference does not exceed 7%. An 
interesting fact is that the agreement between q factors 
obtained by us and by Dzuba and Flambaum [ij] is (on 
average) better than for the states of the Sd^isAp config- 
uration. 



Analyzing the q factors presented in Tables IIIII and IIVI 
we see that the former are (on average) smaller than the 
latter. It has simple explanation because the transition 
between Sd^As^ and Sd^AsAp configurations is basically 
a one-electron 4s-4p transition. The transition between 
3d^As^ and 3(i^4p configurations is due to simultaneous 
4s-4p and As-Sd one-electron transitions. As was shown 
in 23], when a is changing towards its nonrelativistic 
limit a = 0, changes of one-electron energies of s and d 
states are differently directed. It leads to increasing the 
q factors. 

Note that the states nominally belonging to the Sd'^Ap 
configuration are formed, as a rule, as a result of strong 
mixing of Sd^AsAp and 3d^Ap configurations. An admix- 
ture of the Sd^Asip configuration to the 3d'^Ap configura- 
tion leads to opening a strong electric dipole transitions 
from the states of the 3(i^4p configuration to the ground 
state. It allows one to make these transitions observable. 

Using the results obtained for the [Aspdf] , [Bspdf] , and 
[6spdf] basis sets we are able to estimate uncertainties 
of the calculated q factors. We roughly estimate the un- 
certainty as the difference between largest and smallest 
values of the q factors found for the three basis sets men- 
tioned above. Such conservative estimate allows us to say 
that the q factors for the odd-parity energy levels listed 
in Tables Hill and HVl are within 20% accuracy. 



III. ISOTOPE SHIFT CALCULATION IN Fe II 
AND Fe I 

A. Method 

It is known that the energy levels of two isotopes of any 
element are shifted relative to each other. Total isotope 
shift (IS) is usually divided into mass shift and field shift. 
The former is due to nuclear recoil and the latter is caused 
by the finite size of the nuclear charge distribution. For 
light elements the field shift is much smaller than the 
mass shift and we neglect it in our consideration. 

In relativistic approximation the mass-shift operator 
can be represented in the form of expansion in aZ. In 
Ref. [l^l it was obtained the following expression for this 
operator involving first two terms of the expansion over 
aZ 

HmS = JTT X! (P>-Pk - — + («»«»)",] Pk) , (5) 



where AI is the nuclear mass, Pi is the momentum oper- 
ator, rii = Ti/vi, and cxi is the Dirac matrice of the ith 
electron. 

The Eq. ([5]) can be symmetrized over variables of zth 
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and fcth electrons and written in the following form [iJ 



i.k 



X p 



aZ 



(6) 



This equation differs from Eq. ([5]) only by the terms 
~ {aZY . Eq. ([6]) is more convenient for CI calculations 
because it is symmetric in electrons. This allows us to 
add it to the Coulomb integrals. 

Using Eq. ^ we are able to calculate the isotope shift 
in the frequency, 5uj^'^ , of a transition between two iso- 
topes with mass numbers A and A' . Neglecting the field 
shift the expression for 5uj^'^ can be written as 



5uj 



A,A' 



(7) 



NMS 



SMS J 



1 1 



The first term in curly brackets characterizes so-called 
normal mass shift (NMS), while the second is the specific 
mass shift (SMS). In terms of Eq. (O the NMS term is 
determined by the expression with i — k and the SMS 
term is determined by the expression with i ^ k. 

It is worth noting that the use of relativistic MS opera- 
tor is very important for calculating IS for the transitions 
between the FS energy levels of the ground multiplet. An 
account for the relativistic corrections changes drastically 
(up to change of the sign) the values of the IS obtained 
with the nonrelativistic MS operator. 

For transitions from the states of other multiplets to 
the ground state the IS is less sensitive to the relativistic 
corrections. In particular, for such transitions the coeffi- 
cient A:nms can be found with a good accuracy from its 
nonrelativistic expression /cnms = —a;/ 1823. The value 
1823 is the ratio of the atomic mass unit to the electron 
mass. 

Technically the isotope shift of an energy level can be 
found as follows. The operator Hms can be added to the 
many-body Hamiltonian H with an arbitrary coefficient 
A: 



Hx — H + XHms ■ 



(8) 



When the eigenvalue E\ of the Hamiltonian Hx is 
found, the IS correction to the energy can be obtained as 



AE = 



dEx 



d\ 



A=0 



E+x — E- 
2A 



(9) 



The parameter A should be chosen from the conditions of 
numerical stability and smallness of the nonlinear terms. 
In our calculations A/(2M) was put to be 0.001. 



B. Results for Fe I 

We have carried out calculations of the isotope shifts 
for the transitions from the even-parity states of the 



ground multiplet to the ground state and from the odd- 
parity states of the 3(i^4s4p and id'^Ap configurations to 
the ground state using the [4spd/] basis set. In the first 
case for calculation of both fcsMS and /cnms "we used the 
relativistic MS operator given by Eq. ([6]). In the second 
case the specific mass shift fcsMS was calculated in the 
relativistic approximation, while the normal mass shift 
was obtained from the simple formula /cnms = — w/1823 
with use of the experimental frequencies. 

Table |V] presents the results obtained for fcsMS and 
/cnms ■ Analyzing these results we see that the total mass 
shift /cms = fcsMS + ^NMS is small for the transitions from 
the fine-structure components of the ground multiplet to 
the ground state. Note that it is essential to account for 
both terms in Eq. ^ for calculations of these quantities. 
The contribution of the second term in Eq. ^ to fcsMS is 
comparable to the contribution from the main term of the 
mass shift operator. In contrast, relativistic corrections 
to the transition isotope mass shifts of the odd-parity 
states are small and, in principle, can be neglected. 

As it was shown in Ref. [25l| the isotope shifts in Ti ii for 
certain levels are strongly influenced by core- valence cor- 
relations disregarded in our approach. They contribute 
to fcsMS at the level of 50% or even more. Ti ii, as well 
as Fe I, is an element with open d shell. Its main config- 
uration is 3c?^4s. Of course, the core of Ti ii is less rigid 
than that of Fe i and, respectively, the role of core- valence 
correlations should be greater. Nevertheless, we expect 
that a treatment of the core- valence correlations will lead 
to more significant changes of k^MS, than inclusion of the 
relativistic corrections to the interaction between valence 
electrons. 



C. Results for Fe II 

The method of calculation of different properties of 
Fe II is similar to that for the neutral iron. The results 
of calculation of the q factors for a number of the odd- 
parity states were reported in Refs. [13, HI] . All details 
of the method of calculations can be found in [l3| . Here 
we only briefly recapitulate its main features. 

At first we solved the Dirac-Fock equations to find core 
orbitals Is, 3p3/2 and valence orbitals 2td2,/2, 3^5/2, 
4pi/2i 4p3/2- Then we added virtual orbitals, which were 
constructed using the procedure described in section [Til 
The basis set used for seven-electron CI calculations in- 
cluded s, p, d, and / orbitals with principle quantum 
number n < 7 designated as [Ispdf] . Configuration space 
was formed by SD excitations from the configurations 
3d^4s, 3d^4p, and Sd^Asip. 

The method of calculations of the IS remains the same 
as for Fe i. The results of the seven-electron CI calcu- 
lation of the IS in Fe ii are listed in Table |VT1 To the 
best of our knowledge the only calculation of the IS for 
certain transitions in Fe ii was carried out in [l2|. We 
see reasonable agreement between the results obtained in 
this work with the values found in Ref. [131 . 
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TABLE V; Fe I. Experimental and theoretical transition fre- 
quencies Lo of the even- and odd-parity states respective to the 
ground state (in cm~^), fcsMS (in GHz) and feNMS (in GHz) 
are presented. The values are obtained for the [Aspdf] basis 
set. 
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TABLE VI: Fe 11. Experimental and theoretical transition 
frequencies lj of the even- and odd-parity states respective 
to the ground state (in cm~^), fcsMS (in GHz) and /cnms (in 
GHz) are presented, /cms = fcsMS + /cnms. The values are 
obtained for the \7spdf] basis set. 
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and [Gspdf]. Comparing the results obtained for these 
basis sets we could estimate the accuracy of the q factors 
for the odd-parity energy levels listed in Tables IIIII and 
ITVlat the level of 20%. The accuracy of the q factors for 
the transitions from the even-parity states of the ground 
multiplet to the ground state is significantly higher. We 
estimate it at the level of a few percent. 

The certain odd-parity states are astrophysically in- 
teresting because they were observed in quasar absorp- 
tion spectra. Due to strong configuration interaction the 
magnitudes of the q factors vary significantly between 
the states. This makes Fe i a very attractive candidate 
to search for hypothetical a variation. 

We also computed the mass isotope shifts for Fe i and 
Fe II for the [4spd/] basis set. Comparing the results ob- 
tained for Fe II with the values of [I4I we see a reasonable 
agreement between them. To the best of our knowledge, 
the isotope shifts in Fe i were calculated for the first time. 



IV. CONCLUSION 



We have calculated relativistic frequency shifts {q fac- 
tors) for a number of transitions from the even- and odd- 
parity states of Fe i to the ground state. The calculations 
were carried out for the three basis sets [Aspdf], [5spdf], 
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